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Abstract 

Chaotic evolutions exhibit exponential sensitivity to initial conditions. This suggests that even 
very small perturbations resulting from weak coupling of a quantum chaotic environment to the 
position of a system whose state is a non-local superposition will lead to rapid decoherence. How- 
ever, it is also known that quantum counterparts of classically chaotic systems lose exponential 
sensitivity to initial conditions, so this expectation of enhanced decoherence is by no means obvi- 
ous. We analyze decoherence due to a "toy" quantum environment that is analytically solvable, yet 
displays the crucial phenomenon of exponential sensitivity to perturbations. We show that such 
an environment, with a single degree of freedom, can be far more effective at destroying quantum 
coherence than a heat bath with infinitely many degrees of freedom. This also means that the stan- 
dard "quantum Brownian motion" model for a decohering environment may not be as universally 
applicable as it once was conjectured to be. 
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I. INTRODUCTION AND MOTIVATION 



Chaotic classical systems display a phenomenon known as sensitive dependence on initial 
conditions. Two copies of such a system, prepared in nearly identical states (two distinct 
points in phase space, separated by a very small distance), will evolve over time into widely 
separated states. In an idealized case, the distance between the two points in phase space 
grows as e'^*, where A is the largest Lyapunov exponent of the system. This does not happen 
in quantum mechanics. The (oversimplified) reason is that quantum mechanics is linear; 
thus two "nearly identical" states (i.e., states with a large initial overlap) remain nearly 
identical - their overlap is constant under unitary evolution - for all time! However, Peres 
^1 predicted an analogous phenomenon for quantum systems: two nearly identical systems, 
prepared in identical states but obeying slightly different Hamiltonians (i.e., H and H + 6H), 
will evolve into two different states whose inner product decays exponentially in time when 
H is the quantization of a chaotic Hamiltonian. This idea has recently been studied and 
extended by Jalabert and Pastawski Q]. 

Our interest in this behavior stems from a simple model of decoherence, the process 
which most pure states evolve into mixtures due to interaction with an environment y, 

Simple but ge.e.a. ,„odel of decohe.e.ce ,s generated by an ove.a.. Han„lto„ian 
Hs + Hs + Hssj where the interaction Hamiltonian Hg^ is a product of a system operator 
Os and an environment operator O^. When the system is in an eigenstate \sn) of Os 
(with eigenvalue s„), then the environment evolves according to an effective Hamiltonian 
Hs + SnOg. More generally, the system will be in a superposition of the eigenstates of Os, 
in which case the environment experiences a Hamiltonian that is conditional on the state 
of the system. Thus, if the initial (unentangled) state of the supersystem is |\E'5£:(0)) = 
(^^ a„ |s„(0))) |\E'£-(0)) (where |\E'£-(0)) is simply the initial state of the environment), then 
after a time t the state will become [^seif)) = J2n^n\snit)) \en(t)), where \en{t)) is the 
state into which the environment evolves if the system is in state (equation (^). When 
I {en(t)\em{t)) I = 1, no entanglement occurs, the system remains in a pure state, and there 
is no decoherence. When, however, | (en(^)km(^)) I ~ 0, superpositions of |s„) and \sm) are 
transformed into completely decohered mixtures. 

This model is tremendously simplified - it ignores many important phenomena, including 
the effect of the system Hamiltonian. However, it illustrates a key point: the rate at which 
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an environment induces decoherence in a system depends on the decay rate of the overlap 

{t)\em{t)) I , where 



\ej(t)) = exp (^t{H£ + SnOs)?! |^^(0)) 



(1) 



Thus, we expect that the sensitivity to perturbation of a chaotic environment's Hamiltonian 
will result in rapid decoherence even for very weak couplings. This conclusion is encouraged 
by related analytic [3, ^ and numerical 0, 111 . 1 studies. Some of its aspects are also 
beginning to be investigated experimentally 2, ll^ UM- The nature of the chaotic 
evolution seems to be important for this conclusion 



phenomena is still being debated 17, 



I, and the physics of the related 
It is therefore useful to have an exactly solvable 



model that captures some of the features of quantum chaotic evolutions. 

A well-known feature of completely chaotic classical systems is that at every point in phase 
space there exist stable and unstable manifolds - directions along which a cell respectively 
shrinks and grows exponentially. These manifolds fold in phase space, enabling a given region 
to be always stretching in some direction, yet remain within a bounded volume. We examine 
a system that exhibits such an exponential sensitivity to initial conditions, yet is analytically 
solvable: the inverted harmonic oscillator. We note that our model does not exhibit folding; 
we shall comment on the consequences of this shortcoming in due course. We shall also not 
discuss models with mixed phase spaces - our model is clearly too simple-minded for that. 



II. ANALYSIS OF THE UNSTABLE OSCILLATOR 

The Hamiltonian for an inverted oscillator is H = 2m~^^^^Y^' where we have replaced the 
parameter u'^ for a simple harmonic oscillator (SHO) with — A^. Such an "oscillator" does not 
oscillate at all; it has an unstable fixed point aX {x = 0,p = 0}, but accelerates exponentially 
away from the fixed point when perturbed. The price we pay for a solvable system that 
displays such sensitivity is that its phase space is unbounded; the kinetic energy of the 
system grows approximately as E cc e^^*. Although such a system is clearly nonphysical, it 
is an excellent short-time approximation for some real unstable systems. We believe that it 
can reproduce some of the behavior of a chaotic environment. To be specific, our model can 
be thought of as approximating the exponential sensitivity exhibited by the chaotic systems. 
In real chaotic systems, however, in addition to the folding we have already mentioned, local 
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values of chaotic exponents can be quite different from their time-averages. Moreover, the 
directions of stable and unstable manifolds can vary from point to point. Again, our model 
misses this feature. On the other hand, we are encouraged by the fact that this "integrable 
model of chaos" has been used previously in a different context, and that its conclusions 
seem to be confirmed by numerical simulations jo, 

We consider an inverted harmonic environment (IHE) consisting of one such oscillator, 
and couple it to a system consisting of a single SHO with mass m and frequency u. This 
supersystem is linear, so it can be analyzed with the same master equation techniques used to 
treat quantum Brownian motion (QBM) and other linear problems. We will briefly discuss 
the novel method we use to obtain the coefficients of the master equation. 



A. Obtaining the Coefficients of the Master Equation 



We begin with the standard master equation for a linear system coupled linearly to a 
linear environment (as derived in l2f 



221) 



d 



ih 



[x\ p] + ^ \p\ p] + ^ [x, {p, p}] - F{t) [£, p]) 



(2) 



9t -A(t)[x,[x,p]] + /2(t)[x,[p,p]] 

This form is exact when all the terms in if = Hs + + Hss are linear and quadratic 
in positions and momenta, provided that the coupling terms in Hss involve only position 
operators. Here and throughout, we use x and p to denote the position and momentum of 
the system, and Ui and qi for those of the zth environmental degree of freedom. To analyze 
a particular system-environment combination, we need to obtain the specific values of the 
time-dependent coefficients tu^g., 7^fj, F, /i, and /2. Although the IHE has a single degree of 
freedom, we consider an arbitrary linear environment as long as possible in order to make 
the generality of our method explicit. At the end of the general analysis, we will specialize 
to the IHE. 

Since the Heisenberg equations of motion for the system and environment operators 
match exactly the classical equations of motion for the equivalent variables, we can obtain 
the coefficients of the master equation from its classical analogue, the equation of motion for 
the reduced system. This straightforward approach has been examined previously in [2^, 
but primarily in the context of QBM-like systems with infinitely many degrees of freedom. 
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The supersystem is linear and Hamiltonian, so we can write its trajectories as 



m = mm, (3) 

where z is a 2A^-dimensional vector of the form [x,p,yi, qi,y2, ■ ■ ■] and T(t) is a square 
2N X 2N matrix. An equation of motion gives the derivatives of z{t) in terms of z{t) itself. 
We obtain this by first differentiating equation (jS)) to obtain z{t) = Tz{0), then substituting 
the ^(0) obtained by inverting equation Q: ^(0) = T^^z{t). This yields 

m = fT-'z{t) (4) 

This is an equation of motion for the supersystem; it gives the time derivatives of all the 
supersystem coordinates and momenta in terms of their values at time t. However, this is 
not the equation of motion that we need. When we trace over the environment to obtain 
the master equation, we assume that the environment's state at time t is inaccessible; we 
know only its initial state. We need a different equation, one that respects this constraint. 

To obtain an equation of motion that provides x{t) and p{t) in terms of x{t),p{t), and 
the initial state of the environment, we define a new matrix Tp related to T: 

(5) 

I 6ij for i > 2 

As we evolve the quantum state of the system, we presume that at all times we have access 
to knowledge of (1) the reduced density matrix of the system, and (2) the initial state 
of the environment. The corresponding classical state of knowledge is a vector Zp{t) = 
[x(t),p(t), ?/i(0), gi(0), . . .]. This vector can be obtained using 7^: Zp(t) = Tp(t)zp{0). By the 
same process that led to equation (jD), we conclude: 

ipit) = fpTp-'zp{t). (6) 

This yields x{t) and p{t), but in order to obtain the coefficients in the master equation, 
we need the derivatives of higher order powers of x and p. This requires some simplifying 
assumption; throughout this paper, we will choose to assume that all states are Gaussian. 
This is particularly convenient since Gaussian states form a closed set under linear evolution. 
Such states are completely described by linear and quadratic expectation values of x and 
p, so to characterize their evolution we need time derivatives of x^, p^, and xp as well as x 



and p. This is straightforward; we define the symmetric variance tensor Vp = ZpZp, which 
contains all the quadratic combinations of x and p, and transforms as Vp{t) = TpVp{0)Tj' . 
The time derivative of Vp is thus given by: 

Vp = fpT-'Vp + Vp {TpTp-'Y . (7) 

Now, we need to relate these quantities to the coefficients of the master equation. For any 
time-independent quantum operator A, ^ (^AJ = Ti (^pA^. The derivatives of the relevant 
expectation values are obtained from the master equation: 



d (x) (p) 
d{p) 



m 



dt 
d{x^) 

dt 
dt 

d{{x,p}/2) 
dt 



-mul^ (x) - 7,ff ip) + F{t) 

-{{x,p}/2) 
m 

-2muj^ {{x,p}/2) - 2%« + 2F(t) (p) + 2h^ hit) 



(8) 
(9) 
(10) 

(11) 



mu^l {x') + - (P^) - 7.« {{x,p}/2) + Fit) {x) + h'Mt) (12) 



m 



We can apply the preceding classical analysis to the Heisenberg operators, then equate the 
results with those from equations ip iT^ . Equations (jHl lT2|) imply that the matrix TpT~^ 
that gives the derivatives of the Heisenberg operators must take the form 



( 1/m ...\ 

-mJl^ -7,ff Fyx Fqi 



V ; 



(13) 



where the net force Fit) is the expectation value of a force operator: F(t) = 
Fyiif)yi{^) + Fqiit)qiiQ). Thus, F(t) = Fyi (jji) + F^i (qi). Using this matrix in equa- 



6 



tion (|7j) yields the time derivatives of the system variances (x^), {fP') , and {{x,p}/2): 



dt 

d{{x,p}/2) 
dt 



(Vp)n = -(Rp}/2) 
m 

Wu + (Vp)21 



(14) 
(15) 



{x') + ^ iP') - %« ({x,p}/2) + i {{F{t),x}) . (16) 
Comparing this with the results from equations (jHlfT^. we can solve for f\ and /2 as 

flit) = ^ {{{F,p{t)}) - 2 (f) {p{t))) (17) 

Clearly, these coefficients are nonzero only if the force operator F is correlated with (re- 
spectively) p or X. Since F is always expressed in terms of the t = operators of the 
environment (yi(0), gj(0)), these correlations occur because the evolution mixes system and 
environment operators: p{t) = J^i'^^iZiiO) and x(t) = J^i'^u^ii^)- ^y expanding x{t) and 
p{t) in this way and indulging in some tedious algebra, we obtain fi{t) and f2{t)- They are 
most conveniently expressed as the contraction of two tensors, one of which is the initial 
variance tensor of the environment, 



AyfiO) AyMO) 



(19) 



For lack of better notation, we use Ayq = (| {yq + qy)) — {y) (q) throughout. Technically, 
this is a second cumulant, Ayq = {{yq}), but we have adopted this notation because cumulant 
notation is not widely familiar. 

In terms of Vs, the diffusion coefficients are 



N 



1=1 

N 



i=l 



^iFyiT2^2i+l "^j-fgj'^,2i+l 



Fy-T2^2i+2 



FqiT2^2i+2 



^iFyiTi2i+l 1TliFqiTi^2i+l 



2i+2 



FqiTi 



qi-^l,2i+2 



(20) 

(21) 
(22) 
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We now have all the coefficients of the master equation in terms of elements of the two 
matrices Tp and TpT^^^ (note that while equations ()20H2H) are expressed in terms of T 
instead of 7^, they involve only the ffist two rows, which are identical between T and Tp). 
While complete specification of Tp requires complete specification of the systems and their 
couplings, we can simplify the problem further. The underlying physics mandates that the 
matrix 7^ be of a particular form. We first define the 2x2 matrix 

/ 

Mi = (23) 

where rrii is the mass of the ith degree of freedom in the supersystem, is the mass of the 
system, and </'j(t) is a function determined by the form of Hss- We can then write Tp in 
2x2 block form: 

f Mo Ml M2 ■■■\ 



1 



(24) 



1 
1 

V : : : ••/ 

In order to calculate T~^ , we define Dij as the determinant of the 2x2 sub-matrix (7^)[i,2],[jj]- 
Then T~^ can be written simply as 





( 00 




D2i 


D2A 


1 1 


-ms0o 


00 


-Di3 




D12 








D12 
















D12 



\ 



(25) 



By explicit computation, we verify that TpT^ ^ takes the form of equation (fT^ . with the 
coefficients given by: 
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7c: 



^0 



yi 



qi 



^0 - VOVO 
^000 - 0000 



^0 




N I 

< ^ \l m ■ 



i=l 
N 



rrii 



i=l 



rriiFyiCpi miFgi(f)i 



qi Yi 



(26) 

(27) 
(28) 
(29) 

(30) 
(31) 



B. Master Equation for the Inverted Harmonic Oscillator 



We begin with the supersystem Hamiltonian: 



H 



X H y + aJmsTrieXy 

2m. 2 2m, 2 ^ v s e y 



(32) 



The time translation matrix T is obtained by diagonahzing the equations of motion, which 
yields two normal modes. This transformation is characterized by a new harmonic frequency 
uj, a new inverse frequency A, and a mixing angle 6: 



tan 6 



^ ' + A^ - J{n^ + A^f + 4a^ 



2a^ 



(33) 
(34) 
(35) 



Using these quantities, the 7^-matrix can be expressed in the form of equation ((211), where 
the (piit) are given by: 



, , , n ^sinfcjt) . r, ^sinh(At) 

(j)o{t) = cos^e — ^-^ + sm^e- ^ ' 



UJ 



A 



, , , sin26' / sinfcut) sinh(At)\ 



(36) 
(37) 
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Thus, for this case we immediately obtain the parameters of the master equation by direct 
substitution into equations ()26II31|) . All terms share a common denominator, which we 
denote by D: 



cos^ 6 sin^ 6 sin{ujt) sinh(At) + uX (2 cos{ujt) cosh.{Xt) cos^ 6 sin^ 6 + cos^ 6 + sinfj^) 

uj^^ = ^ (^^ cos^ e-X^ sin^ e + ^-^^ [(cu^ - A^) cos(cut) cosh(At) - 2uX sin(cut) sinh(At)]^ (39) 

(a;^ + A^) sin^ 261 / x w, x / x , ^ , . 

%B = ^ [A sm(cjt) cosh(At) - lu cos{ujt) smh(At)] (40) 

Fy = ^ — ^ (cos^ e cosh(At) + sin^ 9 cos{ut)) (41) 

^ If^s i^"^ + A^) sin 26* / 9 ^ -, / , n , 9 ^ / n \ / , n 

F„ = - J — - ^ (cu cos^ ^ sinh At + A sin^ ^ sin cut ) 42 

y nie 2D 

The diffusion coefficients are naturally described in terms of the elements of the tensors 
in equations ()30II31|) . Factoring out their common prefactor /? = sin^ 26^ (cu^ + A^), we 
obtain: 

ff = nieuX/S (cos^ 6 cosh(At) + sin^ 9 cos(cjt)) (A sinh(At) + u sin(cjt)) (43) 

= ujX(3 (cos^ 6 cosh(At) + sin^ 6 cos(cjt)) (cosh(At) - cos(cut)) (44) 

f^y = /3{uj cos^ e sinh(At) + A sin^ 9 sin(cjt)) (A sinh(At) + uj sin(cut)) (45) 

m = A (^ cos^ e sinh(At) + A sin^ 9 sin(cut)) (cosh(At) - cos(cjt)) (46) 
nie 

/l^ = meUjX/3 (cos^ 9 cosh(At) + sin^ 9 cos(cjt)) (cosh(At) - cos(cjt)) (47) 

/f = (3 (cos^ 9 cosh(At) + sin^ 9 cos(cjt)) {00 sinh(At) - A sin(cut)) (48) 

/f = (3{uj cos^ 9 sinh(At) + A sin^ 9 sin(cjt)) (cosh(At) - cos(cjt)) (49) 

/f = A^(cjcos2^sinh(At) + Asin2^sin(tut))(cusinh(At)-Asin(cut)) (50) 
nipUjX 



III. RESULTS AND ANALYSIS 

Having obtained a master equation describing the evolution of a system coupled to an 
IHE (inverted harmonic environment), we now proceed to examine the consequences of that 
evolution. We have two tools for this analysis: on one hand, the master equation and 
its coefficients determine the instantaneous effects of the environment; on the other hand, 
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we can explicitly evolve an initial state of the system to see the time evolution of state 
properties such as entropy and energy. The master equation itself divides naturally into 
two parts corresponding to the two lines of equation the terms in the first line produce 
renormalized unitary evolution (including external damping and forcing terms, which break 
unitarity), while the last two terms are diffusive and responsible for decoherence. We thus 
divide our analysis into three sections, addressing in turn the quasi-unitary portion of the 
master equation, the diffusive terms in the master equation, and the behavior of evolved 
observables. Our primary results are in equations ()63II64|) . where we demonstrate the linear 
growth of entropy in the system at the rate set by the Lyapunov exponent and obtain an 
approximate decoherence timescale that turns out to be logarithmically dependent on the 
coupling. This implies that isolation from chaotic environments is in a sense exponentially 
difficult. In particular, it is even harder to isolate a system from a chaotic environment 
than the many harmonic oscillators of the QBM environment, where the decoherence time 
is approximately quadratic in the coupling strength Thus, the reader who wishes 

to skip straight to the discussion of the implications of quantum chaos for decoherence may 
skim most of sections IIII Al and II1IB| which analyze the master equation in detail. 



A. Unitary evolution 

The first five terms in the right-hand side of equation (j2]), 

^ p] + i p] + T - ^0 ('1) 

are exactly those for the evolution of an isolated harmonic oscillator subject to an external 
force F{t) and a damping force •y^ait). For convenience, although the term %ff[x,{p,p}] 
breaks unitarity, we refer to these terms in the master equation as the "unitary evolution" 
terms. 

Because the coefficients of all the terms except [p^,p] are time-dependent, however, the 
evolution induced by these terms is not necessarily intuitive. The time-dependence of a;^^ 
and 7;,ft (for various values of the parameters) is shown in Figure ^ We see immediately 
that although the coefficients are initially very close to their bare values (cufg ~ cj^, 7cff ~ 0), 
they begin after a certain time to vary dramatically, and finally appear to converge to a 
periodically diverging function of time. Thus, coefficients that we expect to be positive and 
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FIG. 1: Dependence of the master equation coefficients io'^g and on coupling angle 0, system 
frequency uj, and the effective Lyapunov exponent A of the environment. The base configuration 
is cj = A = rUg = rtif. = 1, = 10"'^: plot (a) shows the dependence of o-!l^{t) on 0, plot (b) shows 
the dependence of 7cff(t) on 0, plot (c) shows the dependence of u;^ff(t) on w, and plot (d) shows 
the dependence of uj1^{t) on A. The units of time in all plots are identical but arbitrary. 

well-defined (and which, in the QBM model, settle down after a while to stable values) not 
only take on negative values, but appear at certain times to become infinite. 

The root of this behavior lies in the common denominator of all the coefficients (equation 
(|38|)). This denominator, the determinant of the upper- left 2x2 sub-matrix of 7^, takes 
on both positive and negative values - and thus passes through zero. Since the various 
numerators in equations ()39II50|) do not change sign in sync with the denominator, all the 
coefficients not only vary from positive to negative, but diverge whenever D = Q. While 
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this phenomenon is unexpected and counter-intuitive, we note that it does not indicate 
unphysical behavior. Detailed analysis (omitted here) shows that near a divergence in the 
coefficients, the system state p is forced to assume a form such that the effects of the super- 
operators in the master equation cancel each other out; thus, when the coefficients diverge, 
their effects sum to a perfectly finite value, and the evolution of the state is completely 
physical. 

We conclude from this that the divergences are not symptomatic of any physical phe- 
nomenon, but are rather a consequence of the system-environment paradigm that we have 
imposed on the supersystem. As the system and environment interact with each other, in- 
formation about the initial Xq and Pq of the system is transferred to the environment (and 
vice- versa). At certain times, all information about a particular linear combination of Xq 
and po has been transferred to the environment, in return for which all information about 
some linear combination of and go resides in the system. This is true only for an instant, 
but at that instant p{t) does not uniquely determine p(0), no differential equation for p can 
exist, and the master equation necessarily breaks down. The mathematical reflection of this 
is that the upper left 2x2 sub-matrix of Tp (which specifies the relationship between (x, p) 
and (xo,Po)) is instantaneously non-invertible, because D = 0. 



For t ~ 0, D = 1. Assuming that the coupling is weak (^ ^ 1), we can expand D for 
< t < -2X~^\oge as 



Since the fraction at the end of equation is 0(1), we conclude that D ~ 1 until shortly 
before the ffist divergence occurs; further, the timescale of that divergence is 



While the first divergence may occur later than tc, it cannot occur earlier. Thus, we have 
a natural time scale in the system that divides time into two regions: an initial period, 
when if: < tc', and long times, when t > t^. In the initial period, the unitary coefficients of 
the master equation are well-approximated by their bare values, whereas in the long-time 
regime the coefficients' values periodically diverge and are only indirectly related to their 
bare values. When t ~ tc, the coefficients' values are difficult to characterize. Finally, we 
note that the critical timescale is reflected in the diffusive coefficients as well (see Figure 12]), 



D ~ 



1 + e^e 



(cj^ — A^) sm{ujt) + 2ujX cos{ujt) 



(52) 




(53) 
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but not necessarily in physical quantities obtained from p(t) itself, as we expect from the 
argument that the divergences are not reflected in physical quantities. 

Nonetheless, t^. is a useful time scale to keep in mind because in the initial period we 
can be assured that the renormalized unitary evolution is very similar to the bare unitary 
evolution; thus, the effects of interaction with the environment during this regime will be 
all but unnoticeable on classical scales. Beyond tc, the terms in the master equation that 
govern the evolution of large-scale structures in phase space may depart dramatically from 
their uncoupled values; we cannot be certain. If we take as the criterion for our model's 
relevance to a real chaotic environment that its effects be small on classical scales, then 
relevance is guaranteed for the initial period, but not for long times. Thus, although we will 
examine the long time behavior of the model at times, we regard the initial {t < period 
as a candidate model of decoherence due to a chaotic environment. 

Another perspective on the same conclusion comes from the fact that any physical chaotic 
environment will have a bounded phase space. The operator Og that couples to a system 
operator Os will thus have a bounded spectrum, and the expectation value of the coupling 



spects this constraint, we can consider it a plausible model for that environment; however, 
because the y operator is unbounded and (y^) increases exponentially in time, our model 
will eventually expand into a larger phase space, and the coupling will dominate the overall 
Hamiltonian. The critical time tc indicates roughly when Hse begins to dominate the overall 
Hamiltonian. 

For completeness, we consider briefly the effects of an environment which is not unstable 
- a simple harmonic oscillator or free particle. The simple harmonic environment does not 
have divergences in the master equation coefficients. Instead, uj^^ and 7cff oscillate stably 
around u"^ and (respectively); u^^ is plotted in Figure El The free particle is somewhat 
more interesting; a wave packet spreads slowly (Ax^ oc t) in the absence of a potential, so in 
this sense there is a weak irreversibility. We'll discuss this in more detail when we consider 
entropy, but as Figure El shows, the divergences characteristic of the IHE do plague the 
free-particle environment. However, t^. is no longer logarithmic in the coupling, but rather 
obeys a power-law; thus the free-particle environment can be more effectively isolated from 
the system (like the SHO environment) by making the coupling strength very small. 



Hs£ 



OsOs will also be bounded by some value 




As long as our model re- 
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B. Diffusive terms 



The last two terms in equation 

-Mt)[x,[x,p]] + hit)[x,[p,p]], (54) 

are non-unitary and diffusive, and produce decoherence. If p is transformed to a Wigner 
function W{x,p), then equation (j21) becomes a Fokker-PIanck type equation. Tlie /i or 
"normal diffusion" term is seen to produce diffusion in momentum space according to W (x 
—fi-^W, while the /2 or "anomalous diffusion" term tends to skew the state according to 

oc — /2^^W^ (see [2^, where "anomalous" diffusion was first identified). This is also 
apparent in equations (01121), where /i appears only in while /2 appears only in 

In QBM, the coefficients /i and /2 of these terms equilibrate to constant values or mono- 
tonically decreasing functions after an initial period; in our IHE model these coefficients, 
like those of the unitary terms, vary widely over time. Unlike the other coefficients, /i and 
/2 are dependent on the state of the environment; however, they can each be written as the 
contraction of a 2 x 2 coefficient tensor with the variance tensor of the environment at t = 
(see equations 



fn{t) = Tr 



fyy \fy^ \ I Ay^ Ayq 
fl'^ ] \ Ayq Ag2 



(55) 



2-' n 

The absolute values of the four subcoefficients for each of /i and /2 are plotted in Figure 12 
for the same base parameters as Figure ^ and for three different coupling strengths. Each 
subcoefficient has a prefactor involving the masses of the system and environment, which is 
ignored in these plots; thus, if the system and environment masses are substantially different, 
one subcoefficient may be promoted over another. The most salient feature of these plots 
is that they appear identical - only at times shorter than is any difference visible at 
all between the various coefficients. Examination of equations fj43ll50p confirms that for 
t ^ A^^, all the subcoefficients are proportional to 9'^D~^e^^^. This also explains the sharp 
distinction between short- and long-time behavior in Figure |21 In the short-time regime, 
D ~ 1, and each subcoefficient is well-approximated by / oc e^^*; in the long-time regime, 
however, D ~ cos{ujt + (p), and / oc e'^*sec(ci;t + (p). Thus, on the log-plot of Figure 
13 we see two distinct lines with slopes 2A and A, the second of which is punctuated by 
the periodic divergences seen in all the master equation coefficients. In order to make this 
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Time (au 



FIG. 2: The absolute values of the diffusive master equation coefficients /i and /2 as a function 
of time, plotted on a logarithmic scale in arbitrary units for three different values of the coupling 
angle 0. Since /i and /2 are dependent on the state of the environment, the four tensor components 
of each are plotted here (note, however, that all eight components are nearly identical in this 
case). Plot series (a-d) shows the tensor components of /i, while plot series (e-h) shows the tensor 
components of f2- The base configuration \s oj = \ = = mg = 1. We emphasize that the 
coefficients actually change sign regularly (at each cusp, in fact); in order to use a logarithmic 
scale, the absolute values must be plotted in place of the signed values. 
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FIG. 3: Selected master equation coefficients indicative of the behavior of harmonic-oscillator (plots 
(a) and (c)) and free-particle (plots (b) and (d)) environments are plotted. We plot Jl^j^jJ^ and 
jvy £qj. ijxiaginary A (that is, the environment is a single SHO with frequency — iA) and for A = 0. 



The base configuration is lo 



2' 



A = 4, Q = Y^j, and m<j = = 1. Plot (a) shows iJ^^juP' 



for A = 4z; plot (b) shows Jl^jiJ^ for A = 0; plot (c) shows /f^ for A = 4i; and plot (d) shows 
^yy £qj, ^ _ q_ emphasize that the coupling {B) is much greater than in the cases examined 
previously, yet there are no divergences in iJl^ or 7 for the SHO environment, and the diffusion 
coefficients (represented by remain relatively small and bounded. Divergences do appear in 
the coefficients for the free-particle environment, but their onset time is polynomial in and 
they do not appear in these plots except for the largest value oi Q {t^ /^)- 
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explicit, we examine the short-time (A~^ < t < tc) behavior of /i, which is responsible for 
entropy production (see section ITlI C|) . when the mass ratio menv/^^sys ^ 1- In this limit, /i 
is dominated by /f^; using the approximations sinh(At) ^ 1, 6* <^ 1 and D ~ 1, we obtain: 




The key point in equation ()56p is that the coefficient of the term that produces diffusion 
in momentum (the primary factor in decoherence) increases exponentially with time and A, 
but is only quadratic in the coupling strength. 

The change in the exponent of the / coefficients at t ~ tc is physically relevant as 
well as mathematically sensible. In the short-time regime, the oscillatory dynamics of the 
system and the hyperbolic stretching of the environment proceed largely independently of 
one another; just as the environment induces only minor perturbations in the system, the 
system does not disturb the environment greatly. Thus, the stretching of the environment 
along its unstable manifold is reflected in the system as diffusion in one of the two phase-space 
dimensions. After t^, however, the interaction Hamiltonian begins to dominate the dynamics 
of the system, the overall dynamics is strongly coupled, and the unstable manifold of the 
environment rotates. Diffusion in the system is averaged over stable and unstable directions, 
and so the diffusion coefficients increase only as e'^*. Nonetheless, because diffusion now 
occurs along all directions in phase space, the entropy of p continues to grow at the same 
rate (as will be shown in section UlI C|l . 

C. Entropy of the Reduced Density Matrix 

Having analyzed in detail the IHE master equation, we turn finally to the actual behavior 
of p{t). Because decoherence manifests as entropy production in p, we first calculate the 
Von Neumann entropy S{t) = Tr(p In p) of the reduced density matrix, and examine the 
dependence of S{t) on various parameters. However, entropy production can reflect not 
only the destruction of quantum coherences but also the destruction of large-scale "classical" 
structures in phase space. In order to verify the ability of the IHE to produce decoherence 
without destroying classical structures, we take the expectation value of energy E = (^H^^ 
as a convenient classical quantity, and show that for appropriate initial conditions a dramatic 
rise in 5* occurs while E remains relatively undisturbed. It is worth noting that because 
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all the states we consider have {x) = (p) =0, this is probably an overly strong condition, 
since large (on classical scales) amounts of energy can be added to the system by a simple 
Galilean transformation at t = 0, which changes nothing of the analysis except for adding a 
constant offset to E{t). 



1. Entropy 

The canonical state to examine is a "Schrodinger Cat" state, typically a superposition of 
widely separated coherent states. We have examined the effects of the IHE on such a state, 
but because the state is not itself Gaussian the analysis is quite messy in our ansatz, and so 
we will discuss instead the behavior of a squeezed Gaussian state that is highly extended in 
X. With respect to a diffusive decoherence process, the most significant features of a "cat" 
state are the interference fringes in W{x^p) that lie between the two Gaussian bumps. They 
are extended in x but equivalently narrow in p. A Gaussian state which is highly squeezed in 
momentum has similar small-scale structure, reflected in the strong off-diagonal correlations 
of p. For convenience, we consider such states. 

In general, Tr(plogp) is difficult to compute. For Gaussian states, however, S can be 
easily computed in terms of the state's scaled area in phase space. The Heisenberg uncer- 
tainty relation AxAp > | provides a fundamental unit of phase space area. We thus define 
Oo = |, and define the phase space area occupied by a Gaussian state as 



Ax2Ap2 - (Axp)^ (57) 

and the scaled area as A = a/ao. See III Al after equation ITOl for an explanation of the {Axp) 
notation. Every Gaussian pure state has an area a = oq; mixed states have a > a^. It can 
be shown 2^ that the entropy is given exactly by 

S = -[{A + l) ln(A + 1) - (A - 1) \n{A - 1)] - ln(2). (58) 

A convenient approximation, which is exact for 5 = and is always accurate to within 
1 -ln2 ~ 0.31, is 

5 = ln(A) = iln(A2), (59) 

where the second equality is useful because A? is easily calculable as the determinant of the 
state's variance tensor. This quantity is related to (but not identical to!) the linear entropy 
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^ = 1 — Trp^ because Trp^ = A ^. Thus, the hnear entropy, <; = 1 — A ^, isa first order (in 
A^^) approximation to 5 = In (A). 

2. Analysis of the S{t) Plots 

In Figures EEl S{t) is plotted for a wide range of initial conditions and parameters. The 
initial conditions consist of the squeezing parameter r = ^ and squeezing angle 6 for both 
the system and the environment (large r and 6 = indicates a state extended in position 
and squeezed in momentum, while = ^ implies the reverse). The parameters of the master 
equation include the bare uj and A, the coupling 6, and mass-ratio e = m(./ms. The base 
parameters used in Figures |3]|ni are: = 1, A = 1, = ^, = 1, me = 1, = 4, re = 2, 
6s = 0, and 6^ = 0. Each of the plots varies one of these parameters, except for the fourth 
plot in Figure El in which the entropy for a non-inverted harmonic environment is plotted 
for various values of the coupling angle. 

All but one of the plots demonstrate the same basic behavior: entropy increases linearly 
as S = \t + So, with periodic modulations. The last plot in Figure El shows S(t) for a 
stable environment; the entropy oscillates in time but does not increase irreversibly. Figure 
ini demonstrates the dependence on A, with the intriguing caveat that for A = we obtain 
not constant entropy, but logarithmically increasing entropy. This is explained by the fact 
that for A = the environment is a free particle, whose wave packet spreads as Ay"^ oc t; 
this mild irreversibility produces entropy growth that is only asymptotically zero. For other 
values of A, Ay"^ oc e^'*'*, and entropy grows linearly. 

Varying the squeezing parameter or squeezing angle of the environment (Figure|3)) changes 
only the initial jump in entropy, Sq. Sq is minimal for an unsqueezed environment (re = 1), 
and increases approximately as log(re + r~^). The dependence on 6e is minimal for e = 1; 
for larger or smaller mass ratios, the dependence on the squeezing angle becomes more 
noticeable. It should be noted that the phrase "initial jump in entropy" refers precisely to 
the intercept of S{t) ~ At + in the long time limit; 5*0 can be negative, which means only 
that the linear growth in entropy is postponed for a time —Sq/X. 

Also in Figure EJ it is apparent that the effect of varying r^ and 6s is merely to modify 
the periodic modulation of S(t). Since the system state rotates in phase space (according 
to the dynamics) over time anyway, changing 6s merely changes the phase of the oscillation 
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FIG. 4: Dependence of system entropy on the initial states of the system and the environment. 
The base parameters are u = A = = nie = 1, = = 4, re = 2, 6s = 0, and 9^ = 0. 

Plot (a) shows the dependence of S{t) on r^, the squeezing parameter of the environment; plot (b) 
shows the dependence of S{t) on r^; plot (c) shows the dependence of S{t) on 6^, the squeezing 
angle of the environment; and plot (d) shows the dependence of S{t) on 9^- All entropies are Von 
Neumann entropies: S = Trplnp, and thus are dimensionless. 

in S{t), while determines the shape and amplitude of the oscillation: it becomes more 
dramatic as increases or decreases from 1 (which eliminates the modulation). 

Figures EHSl show the effects of varying uj, A, e, and 6. Both 6 and e affect 5*0 without 
changing the character of S{t) in any other way. This is particularly interesting for the 
coupling 6; regardless of how small the coupling is, the entropy still grows as e''**. In par- 
ticular, the plot for 9 = corresponds to a crossover time of tc > 15, yet entropy begins 
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FIG. 5: Dependence of system entropy on the parameters of the master equation. The base 
parameters for plots (a-c) are uj = X = rUg = rrie = 1, 6* = g^, = 4, rg = 2, 6s = 0, and 6e = 0; 
for plot (d) the parameters are X = Ai, lo = = rrie = 1, = = 4, re = 2, 0s = 0, and 

9e = 0. Plot (a) shows the dependence of S{t) on 6; plot (b) shows the dependence of S{t) on 
ms/rrie; plot (c) shows the dependence of S{t) on uj; and plot (d) shows the dependence of S{t) for 
a stable environment on 9. We emphasize that in plot (d), entropy oscillates but does not increase 
over time (compare with figure inj. The couplings are much larger for plot (d), in order to make S 
perceptible on the same scale as is used for the unstable environment. 

to grow linearly around t ~ 2.5. Changing A, as mentioned earlier, changes the rate of 
entropy production. Finally, the result of varying u is a change in the periodic modulation; 
the frequency of this modulation is of course uj. Of special interest is the curve for uj = 0; 
here, although the periodic modulation has vanished, S{t) rises not as At but as At + log(At). 
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FIG. 6: Dependence of system entropy on the effective Lyapunov exponent. The base parameters 
are lo = rUg = rUe = 1, = = 4, re = 2, Os = 0, and Of. = 0. The curve for A = shows 

logarithmic increase in entropy. 

The explanation is the same as it was for the A = case; when c<j = the system is a free 
particle, and the linear spread of Ax^ contributes a logarithmic term to the entropy. 



3. Conclusions regarding Entropy and Energy 

The linear growth of entropy seen in Figures|3]|niindicates that the IHE model can continue 
to decohere a system to which it is coupled after other models, such as QBM, have ceased 
to produce substantial entropy. However, it does not demonstrate that the IHE can produce 
decoherence faster than QBM. To examine this issue, we first derive a formula for the rate 
of entropy growth. 
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The exact entropy (equation ()58p) is difficult to work with; therefore we approximate 
S{t) with S{t) = ^ InA^ from equation (j^ . Using equations (jHKUl)), we obtain 

= -%,{t)A' + 2 [mAx' - /2(t)(Axp)] . (60) 

Not surprisingly, a positive dissipation coefficient '-f^g causes the state to shrink in phase 
space; however, as we saw in the analysis of the unitary terms, 7cff can be both positive 
and negative, and thus expands the state as often as it shrinks it. We conclude that the 
7cff term contributes only to the periodic modulation of S(t). The effect of the /2 term is 
highly dependent on Axp, which can be positive, negative, or zero. However, because Axp 
oscillates around as the system evolves, this term will contribute primarily to the periodic 
modulation as well. This leaves: 

|(A2) ^ 2mAx' (61) 

Since Ax^ is strictly positive and fi{t) in the short-time regime t < tc is positive (see 
equation (fSUj)). this predicts monotonic growth in A^, and thus S(t). In addition, we can 
use the previous analysis of f\{t) to approximate 

flit) ~ ft2^2g2At (g2) 

(where k depends on u, A, e, and the initial state of the environment, but not on time) and 
integrate equation (|UT|). obtaining 

S{t)) = log {kAx'^)+ log 9 + \t. (63) 

This equation summarizes our most significant results. The amount of entropy produced 
is linear in t and in A, but only logarithmically dependent on the properties of the initial 
state or the coupling strength. We can invert equation (jUHj) to obtain an approximate 
decoherence time t^, which is the time required to produce a certain amount Sd (e.g., 1 bit 
for a Schrodinger Cat state) of entropy from an initial pure state: 

td^j {Sd - log {kAx') - log^) (64) 

The key point is that is inversely linear in A, but only logarithmic in the coupling 6 
or the initial properties of the environment (contained in n). As the coupling becomes 
very weak, then, we expect only moderate increases in the decoherence time. This should 
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be contrasted with QBM, where the decoherence time has a power-law dependence on the 
couphng strength, so that isolating the system from the environment (while still very difficult 



30|) is not quite as hopeless as in the case studied here. 



Finally, we consider the system's energy. If there is no coupling to the environment, then 
E is a constant of the motion; conversely, when E is substantially disrupted, the effects of 
the environment have clearly become noticeable at classical scales. Since the system energy 
is given by 

Eit) = l[rn^u^x^) + ^^{f)y (65) 
we can immediately calculate its time derivative as 

^Eit) = m^ZlSfl (66) 

Clearly E will grow, because of the /i(t) term. However, also contains /i, but multiplied 
by Ax^. Thus, for states that are highly delocalized (Schrodinger cat states), rapid growth 
in S is achieved relative to the growth in E. In addition, E{0) (the initial energy of the 
system) plays no role in thus, if the initial energy of the system is high, the added energy 
E{t) — E{0) will be negligible in comparison for some period of time. We conclude that for 
initial states that have large Ax^ and large E{0) (i.e., superposition states over classical 
length scales), a relatively long period of time exists over which energy growth is negligible, 
while entropy grows rapidly. An example is given in Figure [7| 



IV. CONCLUSION 



We have studied an oversimplified model for a chaotic environment. Our analysis is exact, 
although in the end we make simplifying assumptions in order to get "thumbnail results" 
(such as equation (jHH)). However, the conclusion that we are led to is significant, and we 
expect it to be quite generally applicable; we have shown that an unstable environment with 
a single degree of freedom can produce decoherence more readily than the canonical QBM 
environment, which requires infinitely many degrees of freedom and a much larger Hilbert 
space. Our analysis necessarily examines a system that has an infinite-dimensional Hilbert 
space (since the phase space volume of the IHE is unbounded), but we can imagine a system 
with a bounded phase space that simulates an inverted oscillator arbitrarily well up to a 
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FIG. 7: Comparison of S{t) and E{t) for an extended state. Since E grows as e^^^, we have plotted 
^ ln(£'/ii^o) so that both curves have the same slope asymptotically; thus all quantities plotted are 
dimensionless. The noteworthy regime is t < 5, where the entropy grows steadily but E remains 
virtually constant. The parameters for this calculation are: uj = 10^^, A = 1, rrig = rUe = 1, and 
= gY2- The initial state of the system is squeezed by a factor of 10'^, with the long axis located 
at an angle of ^ to the x-axis; the environment is initially squeezed by a factor of 16 in p. 

certain time. Such a system would produce the same results over the time of interest, but 
would not be analytically soluble beyond that time. 

A brief comment is in order on the principle that enables our "small" (in terms of degrees 
of freedom and available Hilbert space dimension) environment to be so effective. A good 
measure of an environment's decohering efficacy is the amount of entropy it can produce in 
the system, and the rate at which that entropy is produced. For pure initial states of the 
environment, entropy can only be produced through entanglement; the total entropy that 
can be produced is limited by size of the environment's Hilbert space. Thus, in the long 
run, larger environments can decohere more effectively. The rate at which the entropy is 
produced, however, is limited by the rate at which the environment can explore its phase 
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space. Because a collection of harmonic oscillators is a stable system, small perturbations 
due to the state of the coupled system do not induce exploration of a large volume of phase 
space for any one oscillator. The inverted oscillator, on the other hand, can explore its 
volume much more efficiently when perturbed. 

We view this work as the first step in a larger project: understanding the role that chaotic 
systems play in decoherence. These results for the IHE environment clearly show that a 
particular unstable environment can not only decohere a system in a manner markedly 
different from the standard (QBM) models of decoherence, but also display unexpected 
behavior - e.g., the periodic breakdown of the master equation formalism due to singularities 
in the coefficients. While we believe we have sketched accurately the regime in which our toy 
model represents faithfully some aspects of the behavior of actual chaotic environments, the 
clear next step is to examine such environments directly using numerical simulations, and to 
see which of these results for the IHE remain unchanged. Numerical studies carried out to 
date show a .a.ge of behavio. [i HflQ. We hope that the exaetlv solvable .odel 
we have described will aid in the analysis of the relevance of chaos to decoherence. There 
remains another extremely important regime that we do not examine here: environments 
with many chaotic degrees of freedom (the atmosphere of the earth, for instance). 

In addition to these obvious future directions, this work raises related questions about the 
behavior of open systems. The key step in any decoherence process is that of tracing over 
the environment; if this is not done, the state of the supersystem is an entangled pure state, 
not a mixed state. This is commonly justified by the argument that the environment is vast 
- the whole universe, potentially - and this size guarantees that some of the information that 
has been transferred to the environment has been irreversibly lost. Despite this (reasonable) 
justification, explicit models of open quantum systems are usually bipartite; there is a small 
system, and a larger but still small (as compared to the rest of the universe) environment to 
which the system is coupled. If, however, one imagines a set of "concentric" environments 
8^^\ each much larger than the last, as a model for the entire universe, then the environment 
to which the system is coupled {8'^'^^) serves not as an independent environment, but rather 
as a "communication channel" between the system and the greater universe (see 30|,|3l|, also 
0]). In this model, which we intend to investigate, the entropy of the system is not limited 
by the size of the environment (as it is in the bipartite system - environment) model, but 
only by its own size; a small local environment can lead to redundant records of the preferred 
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observables of the system in the rest of the universe. Our analysis of the IHE model indicates 
that chaotic local environments may act as "amplifiers" |32] and, thus, carry information 
away from the system more efficiently than integrable environments. 

Acknowledgments 

The authors would like to thank Diego Dalvit, Juan Pablo Paz and Bill Unruh, for 
discussions and ideas. This work was supported in part by a grant from the NSA. 



[1] A. Peres, Quantum Theory: concepts and methods, vol. 72 Fundamental theories of physics 

(Dordrecht; Boston: Kluwer Academic Publishers, 1995). 
[2] R. A. Jalabert and H. M. Pastawski, Physical Review Letters 86, 2490 (2001). 
[3] W. H. Zurek, Physics Today 44, 36 (1991). 

[4] D. Giulini and et al, Decoherence and the appearance of a classical world in quantum theory 

(Berlin; New York: Springer, 1996). 
[5] J. P. Paz and W. H. Zurek, in Coherent Matter Waves: Les Houches, session LXXII, edited 

by R. Kaiser, C. Westbrook, and F. David (Springer Verlag, Berlin, 2001), pp. 533-614. 
[6] W. H. Zurek, Reviews of Modern Physics 75, 715 (2003). 
[7] W. H. Zurek, Nature 412, 712 (2001). 

[8] Z. P. Karkuszewski, C. Jarzynski, and W. H. Zurek, Physical Review Letters 8917, 405 (2002). 
[9] M. A. Sakagami, H. Kubotani, and T. Okamura, Progress Of Theoretical Physics 95, 703 
(1996). 

[10] H. Kubotani, T. Okamura, and M. A. Sakagami, Physica A 214, 560 (1995). 
[11] P. A. Miller and S. Sarkar, Physical Review E 60, 1542 (1999). 

[12] H. M. Pastawski, P. R. Levstein, G. Usaj, J. Raya, and J. Hirschinger, Physica A 283, 166 
(2000). 

[13] G. Usaj, H. M. Pastawski, and P. R. Levstein, Molecular Physics 95, 1229 (1998). 
[14] P. R. Levstein, G. Usaj, and H. M. Pastawski, Journal Of Chemical Physics 108, 2718 (1998). 
[15] T. Gorin and T. H. Seligman, Journal Of Optics B-Quantum And Semiclassical Optics 4, 
S386 (2002). 

28 



[16] T. Prosen and T. H. Seligman, Journal Of Physics A-Mathematical And General 35, 4707 
(2002). 

[17] P. Jacquod, I. Adagidcli, and C. W. J. Beenakker, Physical Review Letters 8915, 4103 (2002). 

[18] P. Jacquod, P. G. Silvestrov, and C. W. J. Beenakker, Physical Review E 6405, 5203 (2001). 

[19] W. H. Zurek and J. P. Paz, Physical Review Letters 72, 2508 (1994). 

[20] D. Monteoliva and J. P. Paz, Physical Review Letters 85, 3373 (2000). 

[21] P. A. Miller, S. Sarkar, and R. Zarum, Acta Physica Polonica B 29, 3643 (1998). 

[22] S. Habib, K. Shizume, and W. H. Zurek, Physical Review Letters 80, 4361 (1998). 

[23] D. Monteoliva and J. P. Paz, Physical Review E 6405, 6238 (2001). 

[24] P. Bianucci, J. P. Paz, and M. Saraceno, Physical Review E 65, 046226 (2002). 

[25] A. O. Caldeira and A. J. Leggett, Physica A 121, 587 (1983). 

[26] W. G. Unruh and W. H. Zurek, Physical Review D 40, 1071 (1989). 

[27] B. L. Hu, J. P. Paz, and Y. H. Zhang, Physical Review D 45, 2843 (1992). 

[28] J. R. Anglin and S. Habib, Modern Physics Letters All, 2655 (1996). 

[29] W. H. Zurek, S. Habib, and J. P. Paz, Physical Review Letters 70, 1187 (1993). 

[30] W. H. Zurek, in Frontiers of Nonequilibrium Statistical Physics, edited by G. Moore and 

M. Scully (New York: Plenum Press, 1986), vol. 135 of NATO ASI Series B, (quant-ph 

/0212153). 

[31] W. H. Zurek, Physical Review D 26, 1862 (1982). 

[32] R. J. Glauber, in New Techniques and Ideas in Quantum Measurement Theory, edited by 
D. M. Greenberger (New York Academy of Science, 1986), vol. 480 of Ann N.Y. Acad. Sci., 
pp. 336-372. 



29 



